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We apply continuum percolation to proton-proton collisions and look for 

■ the possible threshold to phase transition from confined nuclear matter 
CN ! to quark gluon plasma. Making the assumption that J/tp suppression is a 

good signal to the transition, we discuss this phenomenon for pp collisions, 

■ in the framework of a dual model with strings. 
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■ In recent years, high energy heavy ion collision experiments have been trying 
to collect information on the possible existence of the plasma of quarks and 
gluon (QGP). One of the strategies has been to look for differences in particle 
production between high density central heavy ion collisions, and low density 
ion collisions, nucleon-nucleus collisions and nucleon-nucleon collisions. At the 
CENR/SPS (y/s « 19 GeV) and now at Brookhaven/RHIC (y/s « 130 - 200 
GeV) several important general results have been obtained. 

The charged particle density was found to increases with energy and the num- 
^ ■ ber of participating nucleons. The average transverse momentum also increases 

with energy and particle density, <Pt> increases as well with the mass of the 
produced particle. Strangeness increases with energy and particle density 0. 
All these results, naturally excluding the dependence on the number of par- 
ticipants, are qualitatively similar to results obtained in nucleon-nucleon and 
nucleon-nucleus collisions [[| . None of them, separately, can then be taken as 
clear evidence for the formation of the QGP. 

On the other hand, the anomalous suppression of the ratio J/ip over Drell- 
Yan production ||, at large associated transverse energy Et, has been widely 
accepted as good signal for QGP formation ||. In fact, no such effect was seen 
in lower density nucleus-nucleus, nucleon-nucleus or nucleon-nucleon collisions. 

In this note, we shall argue that if the J/ip suppression has its origin in the 
creation of an extended colour conducting region, as in percolation, the same 
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kind of suppression should occur even in nucleon-nucleon (pp or pp) collisions 
at high enough energy. 

We shall work here in the framework of multi-collision models, namely the Dual 
Parton Model (DPM) || , but try to be as general as possible. The basic ideas 
are: 1) nucleus-nucleus collisions can be built, in a non trivial manner, from 
nucleon-nucleon collisions; 2) nucleon-nucleon collisions occur with formation 
of 2k intermediate strings, strings are always in pairs, k > 1 being a function 
of energy; and 3) strings may fuse and percolate || a process destroying naive 
additivity of elementary collisions. 

The key parameter in transverse plane string percolation is the dimensionless 
tranverse density r], with 



where r is the tranverse radius of the string (we shall take r = 0.2 fm), R 
the radius of the interaction area, and N s the number of strings. Percolation 
occurs, in the R — > oo, N s — > oo limit, for rj > rj c « 1.15. 

As fusion and percolation of strings also occur in nucleon-nucleon collisions it 
is clear that the J/ip over D.Y. ratio will be strongly affected as 77 approaches 
T] c . In the pp(pp) case (|I]) becomes 



where R p is the effective proton radius. As the observed increase of particle 
densities with energy is mostly due to the increase in the number of formed 
strings, k is also an increasing function of energy. Thus 77 may reach rj c and 
become larger than r] c . This implies anomalous J/ip suppression. 

Following previous work |],[7| we treat fusion and percolation of strings as a 
two dimensional continuum percolation problem ||. We performed computer 
simulation by throwing iV discs (of radius 0.2 fm) into a uniform region (of 
radius of order of the radius of the proton, 1 fm) and counting the fraction 
of events f(r]) with percolation. In the R — > 00, iV — > 00 limit f(r]) becomes 
a step function with a sharp change at 77 = r] c . As r/R is not so small, finite 
size effects are important, affecting mostly the slope (a) of the function at 77 c , 
but the value of ?7 C itself. The computer simulation results were fitted by the 
function, 

f( v ) = (l + e-^-^ a Y 1 (3) 
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Fig. 1. Percolation probability as function of the transverse dimensionless den- 
sity r\ for pp collisions geometry. The curve shows the fit with (|3|). We obtain 
a = 0.1666 ± 0.0067 and rj c = 1.3584 ± 0.0116. 

and the following values were found for the parameters: a = 0.1666 ± 0.0067 
and r] c = 1.3584 ± 0.0116 (see Fig. 1) 

Assuming now, as in 0, that the J/?/> production is prevented in a plasma 
of colour charges || and that such situation corresponds to percolation and 
creation of a large conducting area || we obtain for the J/ip over Drell-Yan 
ratio, 



R( V ) = K(l - f( v )) = K 



(4) 



where K w 55 is the value of the ratio at moderate energy. 

The problem now is simply the problem of relating rj to y^, In other 
words, we need to obtain a reasonable estimate for the energy dependence of 
k. This is what we shall attempt now. 

If in nucleus-nucleus collisions <v> is the average number of nucleon-nucleon 
collisions and 2k is the number of strings per nucleon-nucleon collision, the 
average number of strings N s is given by, 



j\r =< u > 2k . 



(5) 



In nucleon-nucleon collisions <u>= 1 and N s = 2k. From ([]]) and ([5]) the 



3 



condition for the percolation transition, in the case of r/R <C 1, is 

77 = r? c = (-^ 2 <z/> 2A; « 1.15. (6) 



By interpreting the NA50 anomalous J/^ suppression at \/s 19 GeV as the 
result of percolation, one can try to estimate the number of formed strings 
in nucleon-nucleon collisions. The basic information is that anomalous sup- 
pression is absent in S-U central collisions, but it is present in Pb-Pb central 
collisions ||. This means, T]s-u < 1> 15 and r]p b _ Pb > 1, 15. 

At relatively low energy, it is known, from NA49 and WA98, SPS experiments, 
that the number of nucleon-nucleon collisions in central nucleus-nucleus colli- 
sions is, roughly, (see, for instance |TT|). 



<„>~1.5-^, (7) 



where N p is the number of participants in a central AB, A<B, collision, 



On the other hand 



R^A^fm. (9) 



By using (0), (|[) and @ in (||) we obtain, from S-U 

fc(Vi^l9)<1.7, (10) 

and, from Pb-Pb, 

k(y/s ^ 19) > 1.6 . (11) 

As we are not so confident with these estimates we shall include an error of 
the order of 15% and study the ratio J/ip over DY in the range 

1.4<fc(v^^l9)<1.9. (12) 

If we look now at the charged particle densities in pp(pp) collisions, in the 
spirit of DPM, we have that particles are emitted from two kinds of strings: 
valence strings(V), always 2 from valence quark interactions and shorter sea 
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Fig. 2. Pseudo-rapidity density as function of cm. energy. The solid line represents 
parameterization used to fit the data. 

strings(S), from sea parton interactions, in a number growing with energy, 
2{k — 1). The central charged particle density is written as(see, for instance, 



TGH and 111 
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dy 
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+ 2(*-l)£ 



(13) 



On the right hand side of fll3|) we have both contributions, from V and S 
strings. 

We assume that ^\v an d ^rls are constant ("Feynman scaling") and that 
the observed rise of the plateau is determined by the increase in the number 
of strings, i.e., by increase of k. In the low energy limit, k — ► 1, we thus have 



Feynman scaling with, from data |TT]JT2 



dN 
dy 



k -> 1 2 



PP 



dN 
dy 



1.45 ±0.05 



(14) 



In order to determine the energy dependence of k we do the following. Solve 



de equation (|I 

terization to the pp(pp) data [11 



for k, with ^\v fixed by 



and using for ^ \ pp a parame- 



dN 
dy 



0.957 + 0.0458 In + 0.0494 In 2 



(15) 
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Fig. 3. The ratio of J/t/j to Drell-Yan events as function of cm. energy in pp collisions. 

For each of the limiting values of k(y/s ~ 19) (|1^) , combined with the limiting 
values of ^\v (|14|), we adjust the constant to obtain agreement with (|i~5|). 
The fit to ^\pp high energy data for y/s > 20 GeV is shown in (Fig. ||). 

In Fig. ^|we present the y/s dependence of the J /if} over DY ratio, equation (f|) 
and (H]), for the two limiting values of k. In conclusion: we expect a fast drop 
of J/i/j over DY ratio in pp(pp) collisions in the energy range 150 < y/s < 1000 
GeV/nucleon. These are energies of RIHC and Tevatron. 

Our work can be criticized from several different points of view: 

1) The J/ip over DY ratio is also affected by internal absorption and as 
the number of strings increases with energy, absorption makes the ratio 
continuously decrease with energy. This effect was not included. 

2) As charm production probability increases with energy the J/ip over DY 
ratio should have a tendency to increase with energy. This correction was 
also not included. One should, perhaps, in future consider the ratio J/ip 
over cc production as the reference quantity. 

3) As in pp(pp) collisions the interaction radius increases slowly with energy 
the critical value of 7], rj c , should decrease with energy. However the im- 
portant effect is the actual decrease of rj ~ with the consequence that 
the percolation transition occurs at higher value of energy. 

4) The model is purely soft model and hard effects related to the increase 
of <pr> and changes in multiplicities were not included. These effects 
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can be accounted for with fusion of strings but were not considered here. 
They are being studied now. 

5) Non-uniform distributions in impact parameter (like gaussian distribu- 
tions) give rise to an increase of r\ c (see last paper in||), and consequently 
the percolation transition will tends to be displaced to higher energy. 

6) One may question the validity of continuum percolation arguments when 
the ratio r/R is so large, r/R^i 1/5. Technically there is no problem, but 
we are not sure of the validity of the treatment. 

While finishing this paper we became aware of the work of T. Alexopoulos et 
al. |13|] dealing with evidence for deconfinement at Tevatron (a/s = 1.8 TeV). 
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